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Abstract 

Among the models of quantum computation, the One-way Quantum Computer [101111] is one of 
the most promising proposals of physical realization [12], and opens new perspectives for paralleliza- 
tion by taking advantage of quantum entanglement [2]. Since a one-way quantum computation is 
based on quantum measurement, which is a fundamentally nondeterministic evolution, a sufficient 
condition of global determinism has been introduced in [4] as the existence of a causal flow in a 
graph that underlies the computation. A 0(n'^)-aIgorithm has been introduced [6] for finding such 
a causal flow when the numbers of output and input vertices in the graph are equal, otherwise no 
polynomial time algorithirQ was known for deciding whether a graph has a causal flow or not. Our 
main contribution is to introduce a 0(n^)- algorithm for finding a causal flow, if any, whatever the 
numbers of input and output vertices are. This answers the open question stated by Danos and 
Kashefi [4] and by de Beaudrap [6]. Moreover, we prove that our algorithm produces an optimal 
flow (flow of minimal depth.) 

Whereas the existence of a causal flow is a sufBcient condition for determinism, it is not a 
necessary condition. A weaker version of the causal flow, called gflow (generalized flow) has been 
introduced in [3j and has been proved to be a necessary and sufficient condition for a family of 
deterministic computations. Moreover the depth of the quantum computation is upper bounded by 
the depth of the gflow. However, the existence of a polynomial time algorithm that finds a gflow has 
been stated as an open question in [3]- In this paper we answer this positively with a polynomial 
time algorithm that outputs an optimal gflow of a given graph and thus finds an optimal correction 
strategy to the nondeterministic evolution due to measurements. 

Keywords: Graph Algorithms, Quantum Computing 

1 Introduction 

A one-way quantum computation [TU] consists in performing a sequence of one-qubit measure- 
ments on an initial entangled quantum state described by a graph and called graph state [8] 
where some vertices correspond to the input qubits of the computation, others to the output 
qubits and the rest of the vertices correspond to auxiliary qubits measured during the compu- 
tation. Since quantum measurements are nondeterministic, a one-way quantum computation 
requires corrections which depend on the results of the measurements, and which should induce 
a minimal depth for the computation. 

Because of these corrections, not all graph states can be used for deterministic computation. 
The measurement calculus [5] is a formal framework for one-way quantum computations, where 
the dependencies between measurements and corrections are precisely identified. Using this 
formalism, Danos and Kashefi in |,4J proved that a one-way quantum computation obtained by 
translation from a quantum circuit is such that the underlying graph satisfies a causal fiow 
condition (see section O) 



*LIG, University of Grenoble, France, mehdi.mhalla@imag.fr 

^Oxford University Computing Laboratory, simon.perdrix@comlab.ox.ac.uk 

^an exponential time algorithm is proposed in [B] 



1 



In [7] a polynomial time algorithm in the size of the graph has been proposed for finding 
a causal flow when the numbers of outputs and inputs are equal, whereas the existence of 
a polynomial time algorithm in the general case, has been stated as an open question. We 
propose in this paper a faster and more general algorithm for finding a causal flow, whenever 
the numbers of inputs and outputs are different. 

It turns out that the existence of a causal flow is not a necessary condition for determinism. 
A weaker flow condition, the gflow condition has indeed been introduced for characterizing 
uniform, stepwise and strong deterministic computation, where the correction strategy does 
not depend on the measurement basis (see [3] for a formal definition.) Here, we introduce 
a polynomial time algorithm for finding a gflow and thus checking whether a graph allows a 
uniform deterministic computation, which gives substantially more relevance to the notion of 
gflow introduced in [3]. 

We also prove that the algorithms proposed are optimal, which means that they give a 
minimal depth flow. This implies that the gflow algorithm gives a lower bound on the complexity 
of a correction strategy in a measurement-based setting for quantum computation. 

2 Definitions 

A graph with input and output vertices is called an open graph, and is defined as follows: 

Definition 1 (Open Graph) An open graph is a triplet {G,I,0), where G = iy,E) is a 
undirected graph, and /, O C y are respectively called input and output vertices. 

During a one-way quantum computation all non output qubits (represented as non output 
vertices in the corresponding open graph) are measured. Since quantum measurements are 
nondeterministic, for each qubit measurement, a corrective strategy consists in acting on some 
unmeasured non input qubits, depending on the classical outcome of the measurement, in 
order to make the computation deterministic. Thus, a corrective strategy induces a sequential 
dependance between measurements. As a consequence, the depth of the quantum computation 
depends on the corrective strategy. 

Corrective strategies will be defined by flows on the open graphs. A flow {g, -<) consists 
in a partial order -< over the vertices [i -< j if i is measured before j) and a function g that 
associates with each vertex, the vertices used for correcting its measurement (all non output 
qubits are measured.) Input qubits cannot be used for correction (see [5].) 

Given an open graph, two kinds of flows are considered: the causal flow and the gflow 
(generalized flow.) The former has been introduced by Danos and Kashefi [4j and corresponds 
to the computation strategy that consists in correcting each qubit measurement by acting on 
a single neighbor of the measured qubit. For a given open graph, a causal flow is characterized 
by a function g which associates with each non output vertex a non input vertex used for the 
correction of its measurement. More formally: 

Definition 2 (causal flow) ((?,-<) is a causal flow of{G,I,0), where g : V{G)\0 V{G)\I 
and -< is a strict partial order over V{G), if and only if 

1. i -< g{i) 

2. if j € N{g{i)) then j = i or i ~< j, where N{v) is the neighborhood of v 

3. i G N{g{i)). 

An example of causal flow is given in Figure 1. Notice that if the numbers of input and 
output vertices are the same, a causal flow can be reduced to a path cover and then to a standard 
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Figure 1: Example of open graph - squared vertices represent inputs, white vertices represent outputs - which 
has a causal flow ((/,-<), where g{ai) — hi, g{hi) = Ci and ao ^ fli ^ 02 -< {&o,bi,62} -< {co,ci,C2}. 



network flow [6]. This reduction has been used to define an 0(n^)-algorithm for finding a causal 
flow in the case where the cardinalities of input and output qubits are the same [6] . 

The second type of flow considered, the generalized flow, gflow, has been introduced in [3] 
and corresponds to a more general correction strategy that associates with each non output 
vertex a set of vertices used for the corresponding correction (instead of a single vertex.) This 
generalization not only leads to a reduction of the computational depth, but also provides a 
corrective strategy to some open graphs having no causal flow. Moreover, notice that gflow 
characterizes uniform, strong and stepwise deterministic computations [3]. 

For a given open graph, a gflow {g, -<) is characterized by a function g which associates with 
each non output vertex, a set of non input vertices used for its correction, and a strict partial 
order -<: 

Definition 3 (gflow) {g, ^) is a gflow of {GJ,0), where g : V{G) \0 ^ p{V{G) {0} 
and -< is a strict partial order over V{G), if and only if 
1- if j e g{i) then i -< j 

2. if j € Odd{g{i)) then j = i or i ^ j 

3. i G Odd{g{i)) 

Where Odd[K) = {n , \N{u) H -fr| = 1 mod 2} is the odd neighborhood of K, i.e. the set of 
vertices having an odd number of neighbors in K. 

A graphical interpretation of the generalised flow is given in Figure HJ 

Past 
• 

Figure 2: Graphical interpretation of a gflow {g, -<): for a given vertex u all the vertices larger than u are in 
the future of u since the corresponding qubits will be measured after the qubit u, all others are in the past of u. 
The set g{u) has to be in the future of u and such that the following parity conditions are satisfyied: there is an 
odd number of edges between g{u) and u and there is a even number of edges between g(u) and any vertex in 
the past of u. 




A flow [g, ■<) of (G, I, O) induces a partition of the vertices of the open graph: 
Definition 4 For a given open graph {G,I,0) and a given flow {g,~<) of {G, 1,0), let 

('max^(F(G)) ifk = 

" \max^(y(G)\y,l,) tfk>0 
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where max^(X) = {u (z X s.t. Vi" G X, -i(n -< v)} is the set of the maximal elements of X. The 
depth d'^ of the flow is the smallest d such that = 0- {yi^)k=o...d-< a partition ofV{G) 
into d^ + 1 layers. 

A causal flow or a gflow (g, -<) of (G, /, O) leads to a corrective strategy for the corresponding 
one-way quantum computation, which consists in measuring the non output qubits of each layer 
in parallel, from the layer to the layer V{^. After the measurement of a layer V^"^, with 
A; > 0, corrections are realised according to the function g by acting on qubits in Uj<fcl^"^ (see 
[3j for details.) The depth of such a one-way quantum computation is d^. 

Definition 5 For a given open graph {G,I,0) and two given flows {g,~<) and {g',~<') of 
{G,I,0), {g,^) is more delayed than {g',^') if\lk, | Uj=o...fc ^fc^l ^ I ^i=o...k I O'f^d there 
exists a k such that the inequality is strict. 

A causal flow (resp. gflow) (g,^) is maximally delayed if there exists no causal flow (resp. 
gflow) of the same open graph that is more delayed. 

For instance, the flow (g, -<) described in Figure 1 is a maximally delayed causal flow. However, 
{g, -<) is not a maximally delayed gflow since {g' , -<') is a more delayed gflow, where ^'(ao) = 
{bo,bi,b2},g'{ai) = ^2}, 5'(a2) = {b2},g'{bo) = {co},g'{bi) = {ci},g'{b2) = {62}, and 
{00,01,02} {^05^17^2} {co,ci,C2}. One can prove that (5', -<') is a maximally delayed 
gflow. 

The following two lemmas are proved for both kinds of flows. 

Lemma 6 If{g,~<) is a maximally delayed causal flow (gflow) of (G, 1,0) then V^^ = O. 

Proof Let {g,~<) be a maximally delayed causal flow (gflow) of {G,I,0). Elements of V^^ 
have no image under g because of condition 1 in both deflnitions thus V^"^ C O. Moreover, by 
contradiction, if O \ / 0, let ^'=^ \{0 \ Vq-^) x V{G). {g, is a causal flow (gflow) of 
{G,I,0): condition 1 of both definitions is satisfied by -<', because the domain of g does not 
intersect O, so for any i in the domain of g, i -<' j iS i -< j; conditions 2 and 3 of both definitions 
are satisfied in a same way. Thus, {g, -<') is a causal flow (gflow) of (G, /, O). Moreover, for any 
k, Uj=o...A:^fc^ ^ ^1=0. ..kVj^ , and \Vf^\ < \V(^ \ thus {g,<') is more delayed than ((?,-<) which 
leads to a contradiction. □ 

Lemma 7 // {g,~<) is a maximally delayed causal flow (gflow) of {G,I,0) then {g,-<) is a 
maximally delayed causal flow (gflow) of {G,I,0 U V{^) where g is the restriction of g to 
V{G) \ (V U ^1^) and ^ =^ \V{^ x y^^. 

Proof First, one can prove that {g,-<) is a causal flow (gflow) of {G,I,0 U V{^). Moreover, 
by contradiction, if there exists a causal flow (gflow) (g' , -<') that is more delayed than (g, ^) 
then it could be extended to (5", -<") where g"{u) = g'{u) if n G F \ {V^f^ U V{^), g"{u) = g{u) 
if li G and -<"=<' U{(n, t>),n G f\u < v}. {g", -<") is then a more delayed causal flow 
(gflow) of {G,I,0) than (g, -<), which leads to a contradiction. □ 

Lemma 8 // (g, -<) is a maximally delayed gflow, then = {n G V \ O, 3K C O, Odd{K) n 
{V\0) = {u}}. 

Proof First, notice that if {g, -<) is a maximally delayed gflow, then for any u G Vi", g{u) C O 
since u ^ v \i v ^ g{u) (condition 1 of deflnition[3l) Furthermore, by deflnition of Vi' , \i u ^ v 
then V G O thus conditions 2 and 3 of deflnition [3] imply that Odd{g{u)) r\ {V \ 0) = {u}. 
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To prove that any ueV\0 such that 3K C O, Odd{K)nV\0 = {u}, u eV{^, we proceed 
by contradiction. We prove that delaying the measurement of a vertex not in satisfying the 
condition permits to create a more delayed gflow. Indeed, let (g, ^) be a maximally delayed 
flow of {G,I,0) and let ui e V \ V be such that 3K C O, Odd{K) nV\0 = {ui}. Let 
g'{u) = K if u = ui and g'{u) = g{u) otherwise. Let ^' be the strict partial order defined by 
u V if u ^ ui and u ^ v or if u = ui and v G K. It leads to a contradiction since {g' , -<') is 
a more delayed gflow of {G, /, O) than (g, -<). □ 

In a similar way, one can prove that: 

Lemma 9 If (g, -<) is a maximally delayed causal flow, then = {u(zV\0,3v(zO, 
N{v)nV\0 = {u}}. 

Lemmas [8] and [9] show that in a maximally delayed flow, all the elements that can be 
corrected at the last step are in the maximal layer of \ O (i.e. in V{^.) Combined with the 
recursive structure of maximally delayed flow (lemma [7]) , this shows that the layers V/^ of a 
maximally delayed flow can be iteratively constructed by finding elements that can be corrected 
starting from the output qubits. This gives rise to the polynomial time algorithms of the next 
sections. 

3 Causal flow algorithm 

The problem of finding a causal flow of a given open graph is presented in [IJ, and a solution 
has been proposed in the case where the numbers of inputs and outputs are the same. The 
complexity of the algorithm is in 0{nm) where n is the number of vertices and m the number 
of edges (more precisely 0{km) where k is the number of inputs (outputs) ^.) We present 
here a more general and faster algorithm. 

Theorem 10 For a given open graph (G, I, O), finding a causal flow can he done in 0{k.n + m) 
operations where n = \V{G)\ is the number of vertices ofG, m = \V{E)\ is the number of edges 
and k = \0\ is the size of the output. 

In order to prove Theorem [TOl we introduce the algorithm [1] which decides whether given an 
open graph has a causal flow, and outputs a maximally delayed causal flow if one exists. This 
recursive algorithm is based on the recursive structure, pointed out in the previous section, of 
the maximally delayed causal flows. 

The algorithm recursively finds the layers (V^"^)fc=o...d^ • &t the k*^ call to Flowaux, the 
algorithm finds the set Vj^ = Out' (see algorithm [1] and figure [3j) To improve the complexity 
of the algorithm, a set C of potential correctors (Vu G Vj^,g{u) G C) is maintained. The 
algorithm produces a subset C' of C of vertices that can be actually used as correctors, the 
set Out' of vertices that can be corrected by C' is produced as well. For the recursive call, the 
vertices of Out' are added to the potential correctors, whereas the vertices used as correctors 
(i.e. G') are removed from the set of potential correctors since a vertex can be used to correct 
at most one other vertex. 

The partial order -< of the flow found by the algorithm is defined via a labeling I which 
associates with each vertex the index of its layer. As a consequence, for any two vertices u and 
V, n ^ f iff l{u) > l{v). 

Proof of Theorem llOt 
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input : An open graph 
output: A causal flow 

Flow (G,I,0) = 
begin 

for all V £ O do 

I /(^;):=0; 
end 

Flowaux (G,I,0,0\I,1); 

end 

Flowaux (G,ln,Out,C,k) = % precondition: Out \ C cannot be used to correct elements of layer k 
begin 

Out':=0; 

C':=0; 

for all V £ C do 

if 3u s.t. N{v) n(V\ Out) = {u} then 
g{u) ■- V ; 
l{v) ■- k; 

Out' ■- Out' U {u}- 
C ■- C U {v}; 

end 

end 

if Out' = then 

if Out = V then 

I true 
else 

I false 
end 

else 

I Flowaux (G,ln,Out U Out',(C \ C')U(Out'nl/\ln),k+l) 
end 

end 

Algorithm 1: Causal flow 



V \ Out := Vi< Out := U«u»-. Vi' 




Figure 3: Causal flow algorithm: At the fc*'' recursive call, the algorithm flnds out the set Vf^ composed 
of the qubits that will be measured at the — k + 1 step of the one-way quantum computation, where is 
the depth of the computation. At that step, all the qubits in Out ~ Ui=o fc-i ^k'' ^^e not measured, whereas 
the qubits in IJi>fc S'l's already measured. The correctors of the elements of V;^ are in a set C C Out of 
candidates composed of vertices not already assigned to the correction of some future measurement. The first 
stage of the algorithm to find out the set V/^ consists in searching, among C, for the elements that have a unique 
neighbor in 1/ \ Out. Let C' C C be this set of correctors. Then, the neighborhood Out' of C' in V \ Out 
is a set of elements that can be corrected at that step, so Out' is nothing but V^!^. For the recursive call, the 
elements of Out' are added to both Out and C, whereas the elements of C' are removed from C. Since at each 
step, a maximum number of vertices are added to V^T', the causal flow, if it exists, produced by this algorithm 
is maximally delayed. 
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By induction on the number of non output qubits, we prove that if the given open graph has 
a causal flow then the algorithm outputs a maximally delayed one. Assume that the given open 
graph has a causal flow. First, if there is no non output qubit, then no correction is needed: the 
empty flow {g, 0) (where 5 is a function with an empty domain) is a maximally delayed gflow. 
Now suppose that there exist some non output vertices, according to lemma [U] the elements of 

satisfy the test at line 13, moreover the precondition at line 8 (that can reformulated as 
9{Vi^) C C) implies that is composed of the elements that satisfy the test at line 13. Thus, 
after the loop (line 19), Out' = and C = g{V{^). Since the existence of a causal flow is 
assumed, cannot be empty (all non output qubits have to be corrected), thus the algorithm 
is called recursively. Lemma [7] ensures the existence of a causal flow in {G,I,0 U V{^) and 
since the vertices in C have no neighbor mV\{OU V{^), they can be removed from the set of 
potential correctors, preserving the precondition. 

The induction hypothesis ensures that the recursive calls output a maximally delayed causal 
flow in {G,I,0 U V{^) and thus the causal flow {g, -<) defined is a maximally delayed causal 
flow of (G,/,0). 

The termination of the algorithm is ensured by the fact that the set of output qubits strictly 
increases at each recursive call. 

For a given open graph, if the algorithm outputs a flow {g,~<), then this flow is a valid 
causal flow since every output qubit has an image under g, moreover for any vertex i, i ^ g{i), 
and finally if j G N{g{i)) then j = i or i ~< j. Thus, if the given open graph has no flow, the 
algorithm returns false. 

To analyze the complexity of the algorithm, we consider the cost for each vertex u, which 
can be decomposed in: 

• Insert u in the set C of potential correctors 

• Create the set of vertices that u might correct N{u) D {V \ Out) 

• Check whether u can be removed \N{u) Ci {V \ Out)\ = 1 

• Update the potential correctors sets N(v) D {V \ Out) for v £ C when u is removed {u 
belongs to Out for the recursive call). 

As there is at most \Out\ potential correctors, the cost for a vertex u can be decomposed 
in: insert in C + create + |OMt| (remove + check). 

Using a data structure for storing the sets N{u) ri{V\ Out) (for example an array with two 
pointers respectively to next and previous elements), one can remove an element in constant 
time and test whether the set contains exactly one element in constant time. For the creation 
of the structure, one needs to compute the intersection of the neighborhood with (V \ Out). 
Checking whether a vertex is in {V \ Out) can be done in constant time by maintaining an 
array of the vertices that are to be corrected, thus given the adjacency list of a vertex u the 
cost of creating N{u) n (^ \ Out) is the degree of u and the total finding cost is 0{m) where 
m is the number of edges of the graph. The overall complexity is then 0{n\Out\ + m). □ 

This result improves the algorithm in [6] that decides, under the precondition |/| = \0\, 
whether an open graph (G,I,0) has a causal flow in 0{km) operations, where k = \0\. In 
[9], Pei and de Beaudrap have proved that an open graph having a causal flow has at most 
(n — l)k — (2) edges. According to this result, the algorithm in [6] can be transformed (see [9]) 
into a 0(/c^n)-algorithm, whereas our algorithm becomes a 0(A;n)-algorithm. 
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4 A polynomial algorithm for gflow 

Theorem 11 There exists a polynomial time algorithm that decides whether a given open graph 
has a gflow, and that outputs a gflow if it exists. 

Proof Let {G,I,0) be an open graph. The algorithm gFlow(r,/, O) (Algorithm [2j) , where 
r is the adjacency matrix of G, finds a maximally delayed gflow and returns true if one exists 
and returns false otherwise. Given a set Out' and a subset X C Out', Ix stands for a |Ont'|- 
dimensional vector defined by Ix(^) = 1 if * £ ^ and Ixii) = otherwise. 



input : An open graph 
output: A generalized flow 

gFlow (r, In, Out) = 
begin 

for all V G Out do 

I ■- ; 
end 

gFlowaux (r, In, Out, Out \ ln,l); 

end 

gFlowaux (r,ln,Out,fc) = 
begin 

Out':= Out \ In ; 
C-0; 

for all ueV\ Out do 

Solve in F2 : rv\0ut,0uflx = ; 

If there is a solution Xq then C:=C U{u} and g{u) :— Xq ; 
l{u) — k; 

end 

if C = then 

if Out=V then 

I true 
else 

I false 
end 

else 

I gFlowaux (r, In, Out U C,k+1) 
end 

end 

Algorithm 2: Generalized flow 

At the k^^ recursive call, the set C found by the algorithm at the end of the loop at line 16 
corresponds to the layer Vf^ of the partition induced by the returned strict partial order. At 
line 13, the columns of the matrix Ty\Out,Out' correspond to the vertices that can be used for 
correction (vertices in Ui<fcy^'^ \ In) and the rows to the candidates for the set Vf^. A solution 
Xq in F2 to Tv\Out,Out'^x = corresponds to a subset of Ui^k^k' \ ^^^^ ^^^^ ^ 
odd neighborhood in Di<kVi^, thus g{u) := Xq satisfies conditions 2 and 3 required by the 
definition of gflow (see Deflnition [3j) . Furthemore, line 10 of the algorithm ensures condition 1, 
thus if the algorithm returns a flow, then it satisfles the deflnition of gflows. 

Now suppose that the graph admits a gflow (g, -<), then it also admits a maximally delayed 
gflow (g' , -<'). The algorithm flnds the set V( (in the loop at line 12), and by induction (similarly 
to the induction in the proof of Theorem llOp it also flnds a maximally delayed gflow in (G, /, OU 
Vi^) with the recursive call. Thus the algorithm flnds a maximally delayed gflow. □ 
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In order to analyse the complexity, notice that lines 10 to 16 consists in solving a system 
Ax = bi for n — i different 6jS where n = \ V\, i = \Out\ and A is a (n — £) x i matrix. In order 
to solve these n — i systems, the {n — i) x n-matrix M = [A\bi . . . bn-i] is transformed into an 
upper triangular form within O(n^) operations using gaussian eliminations for instance, then 
for each bi a back substitution within 0{'n?) operations is used to find Xj, if it exists, such that 
Axi = bi (see pLj). The back substitutions costs O(n^) operations at each call of the function. 
Since there are at most n recursive calls, the overall complexity is O(n^). 

5 Depth Optimality 

We consider in this section the depth of the flows found by the algorithms, which corresponds 
to the number of steps required by the correction strategy, and we show that both algorithms 
find minimal depth flows, hence optimal correction strategies. 

Theorem 12 The previous algorithms find an optimal depth flow 

Proof First, notice that if {g, ^) is more delayed then {g' , -<') then | U-^g Vj^ \ > \ U-^q 
V^'\ = \V\. Thus > d^', = 0, so < d^' . 

Given an open graph (G, I, O), and an optimal depth flow {g, -<), we can define a maximally 
delayed optimal depth flow {g' , -<'): If {g, ^) is maximally delayed then {g' , -<') = {g, -<) other- 
wise {g' , is a maximally delayed flow that is more delayed than (g, -<). According to the pre- 
vious remark (g' , -<') is of optimal depth as well and the optimal depth d{G, I, O) = d^ = d^ . 
By lemma El V{^' = {u, 3K C O, Odd{K) nV\0 = {«}} = V^" where {g" , ^") is the flow 
found by the algorithm. Thus d{G, /, O) = 1 -|- d{G, 1,0 U ). By induction on the number 
of output vertices the algorithm gives an optimal flow for (G, /, O U ) so the depth given 
by the algorithm is the depth of (i^', which is optimal. 

□ 

The optimality of the previous algorithms have several implications in one-way quantum 
computation. First, the depth (optimal or not) of a flow is an upper bound on the depth of the 
corresponding deterministic one-way quantum computation. Moreover, if the one-way quantum 
computation is uniformly, stepwise and strongly deterministic (which mainly means that if the 
measurements are applied with an error in the angle which characterises the measurement, 
then the computation is still deterministic), then the correction strategy must be described by 
a gflow [3j. As a consequence the algorithm 2 produces the optimal correction strategy, and 
the depth of the gflow produced by the algorithm is a lower bound on the depth of a uniformly, 
stepwise and strongly deterministic one-way quantum computation. 

6 Conclusion 

Starting from quantum computational problems (determinism in one-way quanum computa- 
tion), interesting graph problems have arisen like the graph flow dependence of the depth of 
correcting strategies for measurement-based quantum computation. 

We have defined in this paper two algorithms for finding optimal causal flow and gflow. The 
key points are: the simplification of the structure of the flows considering only the maximally 
delayed flows which have a nice recursive structure; a backward analysis (start from the outputs) 
which allows to take advantage of this structure and avoids backtracking. 

From a complexity point of view, an important question is: given a graph state and a flxed 
set of measurements (we relax the uniformity condition) what would be the depth of an optimal 
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correction strategy. One direction to answer this question would be to define a weaker flow 
that is still polynomially computable. 

One can also consider the characterization and the depth of computation in more generalized 
measurement-based models where other planes of measurements are allowed. 
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